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Abstract. Recent interest in discrete, classical integrable systems has focused 
on their connection to quantum integrable systems via the Bethe equations. 
In this note, solutions to the rational nested Bethe ansatz (RNBA) equations 
are constructed using the "completed Calogero-Moser phase space" of matrices 
which satisfy a finite dimensional analogue of the canonical commutation rela- 
tionship. A key feature is the fact that the RNBA equations are derived only 
from this commutation relationship and some elementary linear algebra. The 
solutions constructed in this way inherit continuous and discrete symmetries 
from the CM phase space. 



1. Introduction 

Let M. n be the set of all pairs of n x n complex matrices (X, Z) which satisfy 
the commutation relationship 

(1) rank([X,Z] + 1) = 1 

where / is the n x n identity matrix. (In other words, this says that [Z, X] w /, 
up to a rank one deformation.) Then M. n forms the "adelic Grassmannian" or 
"completed Calogero-Moser phase space" that arises in the study of integrable 
systems of mathematical physics 13 |9|, [ll] as well as in the study of properties 
of rings of differential operators jyT Hj. It is the purpose of this note to establish a 
relationship between these matrices and the integrable system 

2 TT k / 3 --^ -t-t— = -1 Vl< j < n 

for functions X™ (1 < j < n) of the discrete "time" parameter m. Note that equa- 
tions (^|) arise independently as the Bethe equations for the elementary energies of 
certain solvable quantum models and as a discretization of known classical particle 
systems §01. 

One philosophy behind the recent interest in discrete integrable systems is that 
things seem to be simpler, and theorems true for more elementary reasons, in 
these situations. This can be seen for example from the fact that the connection 
between the commutation relationship (|l|) and the system (||) is proved by only a 
few elementary identities of linear algebra. 

2. Some basic identities in the case det X = 

Fix a choice of (X, Z) £ A4 n such that det(X) — and let e and / be vectors such 
that e • f 1 is the rank one matrix [X, Z] + 1. We will use the notation M to denote 
the matrix of cofactors of the matrix M. (In particular, if M is invertible then 
M = det(M)Af _1 .) Any matrix with determinant zero has a matrix of cofactors 
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with rank at most one, so once again we may specify that v and w are two vectors 
such that X = v ■ w T . 

It will be convenient for us to refer to the following functions 

p(X) := P ■ (X~Z) ■ v 

q(X) := w T ■ (XI^Z) ■ a 
and the constants 7 := (w T ■ e) and /i := (Z 1 " • v). 

Lemma 2.1. We have: 

1. det[(AJ -Z)-X + I}= 7p(A), 

2. det[Jf ■ (A/ - Z) - 7] = -jug(A), 

3. det[(AiJ -Z)-X- (X 2 I -Z) + (A 2 - Ai)J] = (A 2 - X 1 ) P (X 1 )q(X 2 ). 
Proof. By the commutation relationship one has that 

(XI - Z) ■ X + I = X ■ (XI - Z) + e • Z 1 ". 

Then, it is a general fact that the determinant of a rank one perturbation of a 
matrix is given by the formula 

det(M + e • f 1 ") = dct(M) + f ■ M ■ e. 

Applying this fact gives the first identity above. The others follow from the same 
sort of argument. □ 

3. Nested Bethe Ansatz Equations 
For any pair (X, Z) e M. n define the function 

00 

Tx,z&X)=T(lX) :=dct(X + J2e i (X l I-Z)- 1 ) 

i—l 

where I = (£i,£2, ■ ■ ■) and A = (Ai, A2, • ■ ■ ) A, S C, only a finite number of 
£j 7^ 0). This is a polynomial r- function for the KP hierarchy written in so-called 
Miwa variables. One could show that in any three of the variables li this function 
satisfies the discrete Hirota equations Q. However, here we will only observe that 
the linear algebra identities above show that the roots of r satisfy the RNBA 
equations (§). 

The RNBA equations involve only two variables: x and m. These can be taken 
to correspond to any two pairs (£i,Xi). So, let us choose any two values for i: i\ 
and i2- We adopt the notation of Q, specifying a non-zero lattice spacing 77 G C, 
and writing the function r as r m (x) where x — r)li x and m = £i 2 G C. 

Note that T m (x) is a polynomial in x (a special case of the elliptic polynomials 
considered in (Jt| ) . Denote by xj 1 the roots of the polynomials r m (x) (listed in any 
order). The connection between these functions and equations (|^) is an immediate 
consequence of the following factorization of r evaluated at specific values of its 
arguments: 

Lemma 3.1. Choosing either the top or bottom sign at every choice, we have the 
two factorization formulas: 

r m ^(xf ± r,) = —(±x V2 T K) T m ^(x^)T m (xj ± I,). 
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Proof. Let us denote by X' the matrix 

x m 

X' = X + -^-(XiJ - Z)- 1 + m(A 42 - Z)-\ 
V 

Then, since r m (x™) = we know that det X' = 0. Also, it is clear that rank([X', Z]+ 
I) = 1. Hence we may apply the results of the lemma of the previous section to 
the function 



r^ixf + v) = dct(x' + (x h i - zy 1 - (\ t2 i - zy 1 ) 

i 



det(A il J - Z){X l2 I - Z) 



det((A n J -Z)-X' ■ (\ l2 -Z)- (\ l2 - X h )I) 



rP(Ai)g(A 2 ) 



det(X h I - Z)(X i2 I - zy 

Note that we first used part 3 of the previous lemma to factor r m_1 (a;J l + n) into 
a constant times the polynomials p(Xi) and (/(A2) and then used parts 1 and 2 to 
view those as r evaluated at other points. The other equation in which the lower 
sign is chosen at each place is proved similarly. □ 



Corollary 3.1. The functions xj 1 satisfy the RNBA (Q) equations. 

Proof. Using the last lemma, we can factor one term and recombine two other terms 
in this three term product to deduce that 

T m+l (xJ)T m (xJ - ri)T m -\x? + r?) 



= T-'^[xJ)r m {xJ-r,) ( —(X n -X^r^ix^ixj+rj) 



Xi~ ~ A 



7/! 



^ m+1 ( X f)T m {xf - Tj) ) T m - 1 (x j )T m { Xj + V) 



= --- ■ L {xJ - ii)7^ 1 (x J )r'"(i j + ?/). 



So, in particular, 

T m+1 (xT)T m (xT - vW^HxT + n) 



T m+1 (x™)T m {xf - ^t" 1 - 1 ^™ + 77) 
T m+1 {xf - ri)T m ~ 1 {x. ] )T m {x j + 77) 



Expanding this equation by the formula r m (x) — cY\{x — xj 1 ) (for some constant 
c) yields exactly the RNBA equations. □ 

As a consequence, we get our main result: 

Theorem 3.1. Let (X,Z) £ A4 n , then for any Ai,A2 G C (X2 not an eigenvalue 
of Z ), the n eigenvalues of the matrix 

X{m) := -rjX ■ (Ai - Z) - mr;(A 2 - Zy 1 • (A a - Z) 

satisfy equations @). Note that X(m + 1) — X(jn) is a constant matrix and hence 
we have represented solutions of (|2|) as a free flow on matrices. 
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3.1. Remarks. As a result of this construction, several symmetries of the solution 
space manifest themselves. In particular, in addition to the obvious symmetries 
presented by translation in I (the "commuting flows" which appear as translations 
X — > X + f(Z) on M n ), there is also the "dual" flow (cf. |^|, |j) corresponding 
to translations of the form Z — > Z + f(X). In addition, the parameters A allow 
for continuous deformation of the solutions and the involution (X, Z) — > (X T , Z T ) 
on M. n provides discrete symmetry of the of solutions. It would be of interest to 
observe whether this symmetry has any special significance for this system. Recall 
(cf. ^ m |9[ [IT]) that this same involution is a linearizing map for Calogero- 
Moser type particle systems and a "bispectral involution" for the Lax operators of 
rational solutions to the KP hierarchy JL0[ . 

As we know a great deal about the structure of A4 n (jll), this may be useful 
in finding new features of the Bethe ansatz equations. For instance, note that the 
matrix X(m) corresponding to the simplest 2x2 matrices satisfying ([!]): 



X 



on „ _ o o 
oo z - 1 o 



has eigenvalues 



with 



—7] (Aim + A 2 (-1 + Ai + m)) ± er(Ai, A 2 ,to) 
2Xi 



cr(Ai, A 2 , m) = (-4AiA 2 m(-l + A 2 + to) + (Aim + A 2 (-1 + Ai + to)) 2 ) 

So, in particular, in the case Ai = A 2 , they become simply 

{— ?7(Ai - 1 + to), -rjm}, 

a sort of "bound state" in which the eigenvalues remain a distance of rj(Xi — 1) 
units apart regardless of the value of the discrete time parameter to. 

That the manifold M. n should be related to solutions of the RNBA equations 
is not a surprise. In fact, one knows from previous results that these matrices 
can be used to write r-functions, that when written in Miwa form would satisfy 
the Hirota bilinear difference equation closely related to the RNBA equations (cf. 
|6) and [jy]). The point of the observations above is merely that one is able to 
precisely determine the connection between the commutation relationship ([!]) and 
the dynamical system (||) without refering to any of these more general results. 
In particular, the results above provide a means for demonstrating the role of the 
manifold M n in integrable systems using only elementary linear algebra, without 
any mention of symplectic geometry or Hamiltonian dynamics (as was used in jsj). 
It might be of interest to determine to what extent one may rederive the known 
results about Calogero-Moser systems in the continuum limit using only the linear 
algebra identities recalled here. 
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